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Abstract 

This is the first Monte Carlo study of the hard-sphere lattice gas with near- 
est neighbour exclusion on the body-centred cubic lattice. We estimate the 
critical activity to be 0.7223 ± 0.0003. This result confirms that there is a 
re-entrant phase transition of an antiferromagnetic Ising model in an external 
field and a Blume-Emery-Griffiths model on the body-centred cubic lattice. 



1 Introduction 



For phase transitions and critical phenomena of a hard-sphere lattice gas [[T], 
|^. H] while many authors have investigated two-dimensional systems |IJ H, 
H [7|, |8|, |9|, [lO], [11], 0, [L3|, [14], [L5], [Tj|, there are only a few studies in three 



dimensions || [L7], |18[ . One of the reason is that series expansions or transfer 
matrix methods have been applied to the systems mainly. We have obtained 
the critical activity and the critical exponents of the hard-sphere lattice gas 
on the simple cubic lattice by using a Monte Carlo method and a finite-size 
scaling [|17], [18|]. In this paper we perform Monte Carlo simulations on the 
body-centred cubic lattice for the first time and estimate the critical activity. 

Consider a hard-sphere lattice gas whose atoms occupy sites of a lattice 
and interact with infinite repulsion on a site and of nearest neighbour pairs. 
The grand partition function is 

S v (z)=J2z N Z v (N), (1) 

N 

where z is an activity and Zy(N) is the number of configurations in which 
there are N atoms in the lattice of V sites. At z = +oo a configuration of the 
ground state is that the atoms occupy all the sites of one sublattice and the 
other is vacant. There is no atom at z = 0. A continuous phase transition 
occurs at a critical activity. 

The hard-sphere lattice gas relates to an antiferromagnetic Ising model 
in an external field [0, |S], |TT| . The critical curve of the latter behaves as 

H = H c + a*k B T 

near T = 0, where T is the temperature; H is the external field; k-g is 
Boltzmann's constant; H c is the critical field at T = 0. In the present 
paper we shall consider only the case H c > 0. When a* is positive, a re- 
entrant phase transition occurs. The system is in the paramagnetic, the 
antiferromagnetic ordered, and the paramagnetic phase as T is decreased 
when H is fixed slightly above H c . The slope, a*, is given by the critical 
activity, z c , 

a' = — in z c . 
2 

On the body-centred cubic lattice Gaunt || estimated z c to be 0.77(5) by 
using series expansions of the hard-sphere lattice gas and Landau |L9[ ob- 
tained a* = 0.160(16), i.e., z c = 0.726(23) from Monte Carlo simulations of 
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the antiferromagnetic Ising model. These results consist each other within 
errors. There is a re-entrant phase transition since a* is positive. 

In the next section we define physical quantities measured. In section |3] 
we present Monte Carlo results. A summary is given in section |j. 



2 Monte Carlo simulations 



We use the Metropolis Monte Carlo technique |20], |2l| to simulate the hard- 
sphere lattice gas (JJ) on the body-centred cubic lattice of V sites, where 
V = 2xLxLxL (L = 2 x n, n = 2, 3, ... , 12), under fully periodic 
boundary conditions. A body-centred cubic lattice is made up two Lx L x L 
simple cubic lattices. According to Meirovitch |fT3 |, we adopt the grand 
canonical ensemble. The algorithm is described in the references [|T3, |T7| ]. 

We start each simulation from a large activity and then gradually decrease 
an activity. The initial configurations have been obtained from preliminary 
simulations. The pseudorandom numbers are generated by the Tausworthe 
method ^2, 23fl . We measure physical quantities over 10 6 Monte Carlo steps 
per site after discarding 5 x 10 4 Monte Carlo steps per site to attain equi- 
librium. We have checked that simulations from the ground state configu- 
ration and no atom one gave consistent results. Each run is divided into 
ten blocks. Let us the average of a physical quantity, O, in each block (O)i] 
i = 1, 2, . . . , 10. The expectation value is 

i io 



(O) = T7) 5>>*. 
1U i=l 



The standard deviation is 



a(o) = ({oy - {off 2 /V9. 



Let us define a density by 



p = N/V 

where N is the number of the atoms in the lattice of V sites and an order 
parameter by 

R = 2 (N A - N B )/V 
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where (Nb) is the number of the atoms in the A (B)-sublattice and 
N = N\ + A^b- We measure the isothermal compressibility: 



K = V{(?)-(p)*)/<P), (2) 
the staggered compressibility: 

X t = V {(&) - (|i2|> 2 )/4, (3) 

and the fourth-order cumulant of R pl |: 



U = \- (4) 

3 Monte Carlo results 

Figure 1 shows the activity dependence of the isothermal compressibility, 
Kl(z), defined by @ for various lattice sizes. The solid curves are obtained 



by the smoothing procedure of the .B-spline ||25|| . As L increases, the shape 
of the curve becomes sharper. There are shifts of the peak positions. In 
figure 2 we show the activitiy dependence of the staggered compressibility, 
Xl{ z )i defined by (||) for various lattice sizes. It does not seem that the 
position of the peak shifts to contrast those of Kl(z). We show the activity 
dependence of the fourth-order cumulant, Ul(z), of R defined by (f|) for 
various lattice sizes in figure 3. There is an intersection between the curves 
with the size L and L + 2. The positions of these intersections are within a 
narrow region. 

We define effective critical activities, z^ &x (L) and ^ ax (L), as the peak 
position of k l (z) and Xl( z )> respectively, and ^ oss (L) by 

U L {zl oss {L)) = U L+2 {zL ss {L)). (5) 

They will converge to the critical activity, z c , as L — > +oo. We plot them 
against 1/L in figure 4. 

We decide to estimate z c from z^JJu) by the following reasons. Although 
z^ ax (L) seems to behave systematically for L > 6, it is difficult to extrapolate 
z c from it since we need a precise value of a critical exponent v: z^ ax (L)—z c ~ 
L~ x l v . We do not know the value of this system. We cannot see systematic 
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behaviour in z^ ax (L) and z^! TOSS (L) for L > 8. In the latter the value of 
z^ TOBB (20) is deviate from the others with L > 12. We can see from the 
definition (|Sp that it relates to the values of z^ TOSS (18) and z^ TOSB (22). Thus 
we do not adopt z^ mss (L) as an estimator of z c since we cannot ignore it 
simply. We get the result, z c = 0.7223(3), by the arithmetic mean from 
the data z^ ax (L) with L — 12, 14, . . . , 24. It consists with previous results, 
z c = 0.77(5) H and z c = 0.726(23) [|TT| [19| , within errors. Our result is more 
precise than theirs. 



4 Summary 

We performed the Monte Carlo simulations of the hard-sphere lattice gas 
with nearest neighbour exclusion on the body-centred cubic lattice under 
fully periodic boundary conditions. We estimated the critical activity, z c , to 
be 0.7223(3). It consists with z c = 0.77(5) by Gaunt §. 

As is described in section [1] the critical activity relates to the slope of 
the critical curve of the antiferromagnetic Ising model in the external field. 



Our result, a* = 0.1627(2), agrees with a* = 0.160(16) by Landau |]|. The 
system exhibits a re-entrant phase transition since a* is positive. 

In closing this paper we want to mention that an antiferromagnetic Ising 



model in an external field is equivalent to a Blume- Emery-Griffiths model |26 
Kasono and Ono |27| confirms that there is a re-entrant phase transition of 
the latter on the body-centred cubic lattice. Our result supports theirs with 
high precision. 
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Figure captions 



Figure 1 Activity dependence of the isothermal compressibility, k, defined 
by (H) of the hard-sphere lattice gas fll]) on the body-centred cubic 
lattice of V sites under fully periodic boundary conditions. V = 2 x 
Lx Lx L; L = 4: O, 6: +,8: □, 10: x, 12: A, 14: *, 16: o(small), 18: 
• (small), 20: o (middle), 22: •(middle), 24: o (large). The solid curves 
are obtained by the smoothing procedure of the fourth-order 5-spline. 
The data with L = 4, 6, . . . , 14 have been omitted to preserve the 
clarity of the figure. 

Figure 2 Activity dependence of the staggered compressibility, x\ defined 
by (|3p. The meaning of the symbols and the curves is the same as in 
figure 1. The data with L — 4, 6, . . . , 20 have been omitted to preserve 
the clarity of the figure. 

Figure 3 Activitiy dependence of the fourth-order cumulant, U, defined 
by (f|). The meaning of the symbols and the curves is the same as 
in figure 1. The data with L — 4, 6, . . . , 20 have been omitted to 
preserve the clarity of the figure. 

Figure 4 Size dependence of the effective critical activities, z^ ax (ij): +, 
4l x (i): °, and z^ ross (L): x. Errors are less than the symbol size for 
^max(-^) an d ^ctoss(^)- The horizontal line denotes z = 0.7223. 
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